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ABSTRACT 

Sakai et al. have observed long-chain unsaturated hydrocarbons and 
cyanopolyynes in the low-mass star-forming region L1527, and have attributed 
this result to a gas-phase ion-molecule chemistry, termed "Warm Carbon Chain 
Chemistry", which occurs during and after the evaporation of methane from 
warming grains. The source L1527 is an envelope surrounding a Class 0/1 pro- 
tostar with regions that possess a slightly elevated temperature of ~ 30 K. The 
molecules detected by Sakai et al. are typically associated only with dark molecu- 
lar clouds, and not with the more evolved hot corino phase. In order to determine 
if L1527 is chemically distinct from a dark cloud, we compute models including 
various degrees of heating. The results indicate that the composition of L1527 
is somewhat more likely to be due to "Warm Carbon Chain Chemistry" than 
to be a remnant of a colder phase. If so, the molecular products provide a sig- 
nature of a previously uncharacterized early phase of low mass star formation, 
which can be characterized as a "lukewarm" corino. We also include predictions 
for other molecular species that might be observed toward candidate lukewarm 
corino sources. Although our calculations show that unsaturated hydrocarbons 
and cyanopolyynes can be produced in the gas phase as the grains warm up to 
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30 K, they also show that such species do not disappear rapidly from the gas as 
the temperature reaches 200 K, implying that such species might be detected in 
hot corinos and hot cores. 

Subject headings: ISM: molecules, ISM: individual (L1527), stars:formation 



Introduction 



Recent observations of the warm gas associated with L1527 by lSakai et al.l (j2008a[ ) (hereafter 
SSHY) have revealed an unexpected chemical composition including unsaturated chain rad- 
icals of formula C n H, carbenes of formula C n H 2 , and cyanopolyynes, HC2 n C N. These 
molecules are typically associated with the cold, dark cloud stage of star formation J Suzuki et al. 



1992 



Markwick. Millar, fc Charnley 2000 ; Hirota. Maezawa. fc Yamamoto!l2004 ; Smith. Herbst. &: Chang 



20041 ). but independent observations indicate that the physical state of L1527 is more dense 
and slightly warmer than dark clouds. Indeed, the column densities of 13 species detected 
by SSHY toward L1527 are only slightly lower than those detected toward the dark cloud 
TMC-1. The presence of such unsaturated species indicates that L1527 may be one example 
of a chemically unexplored, transitional phase of low-mass star form ation. The analogy with 
dark clouds is streng thened by the detection of the anions CrH~ (ISakai et al.l 120071 ) and, 
most recently, C 4 H" (ISakai et al.ll2008bl : kgundez et alihoosh . 



The physical structure of L1527 contains an envelope surrounding the low- mass pro- 
tostar IRAS 04368+2557 with the object between the Cl ass and I protost e llar p hases 
( Andre. Ward-Thompson fc Barsonyl boool ). According to Myers fc Mardonesl (1998 ). the 
envelope is undergoing inside-out collapse in good agreement with the IShul (119771) model. 



The source also shows evidence for bipolar outflows in the plane of the sky (jOhashi et al. 



19971 ). The SSHY observations include a description of a profile mapping technique, in- 
dicating that the velocity of C4H is consistent with the central protostar rather than the 
outflow region. This region differs from previously investigated star- format ion phases be- 
cause the envelope is only sli ghtly elevated in temperature . Based on a fit to central peak 
emission and the SEP profile, iHogerheijde fc Sandelll (120001 ) determined a temperature of 35 
K, while lLoinard et al.l ( 20021 ) obtained a temperature of 22 K via a similar computation for 
the central region of LI 527. These temperatures are slightly warmer than a typical 10 K 
cold cloud but colder than a hot core/corino region. Selected temperatures in these evolved 
regions have been estimated as T = 50 — 80 K for the low mass hot corinos in NGC1333 



(IMaret et al.ll2004h . T > 100 K for the IRAS 16293 hot c orino (ICeccarelli et al.ll2000h . and 
T = 150 — 200 K for the G327.3-0.6 high mass hot core JGibb et al.lboooh . In addition to 
observations of its physical state, L1527 was included among 18 dense cores in the molecu- 



- 3- 



lar survey of I J0rgensen. Schoier. fc van Dishoeckl (120041 ) (hereafter JSV04). These authors 



detected the molecular species CO, CN, CS, HCN, HNC, HC 3 N, SO, HCO+, and N 2 H+ 
towards L1527. Although this set of observations does not include hydrocarbons or complex 
cyanopolyynes and does not include velocity variations, it does provide an independent basis 
for model comparison. It is clear that further interferometric study of L1527 will provide 
much more information toward understanding the nature of this source. 

While the presence of non-saturated hydrocarbons is more typically a dark cloud char- 
acteristic, the hot core/corino evolutionary phases are more commonly associated with more 
saturated, oxygen- and nitrogen-rich species. What can explain the dark-cloud species de- 
tected in L1527? One possibility is that these molecules are merely remnants of the colder 
pre-stellar core phase that survive a very fast collapse towards the center of the object. In 
the dense portions of such sources, however, it is expected that the carbon chain molecules 
would freeze out to grain surfaces more rapidly than in a cloud of lower density. As noted in 
SSHY, the slight elevation in temperature is capable of regenerating abundances of carbon 
chain molecules after they have frozen out to grain surfaces. The mechanism for this pro- 
cess is more likely the selective release and subsequent reaction of highly volatile molecules 
such as methane (CH 4 ) from grain surfaces than the actual sublimation of carbon chains 
formed by cold processes. For this reason, the process was named 'Warm Carbon Chain 
Chemistry", or WCCC (SSHY). 

In this paper, we present detailed models of a low-mass star-forming region with a 
slightly elevated temperature, for which we use the phrase "luke-warm corino", where the 
adjective "luke-warm" was previously used in talks by E. van Dishoeck. Specifically, we 
attempt to determine if the observations can be explained by models of a cold region alone, 
or if the chemistry is the signature of a newly discovered protostellar phase that requires 
only a small amount of heating. In particular, we consider models in which a cold dense 
phase is followed by a warming to temperatures of 30 K and beyond, based both on the 
dust warming an d on molecular lin e obse rvations (SSHY). The method used is similar to 



that discussed by iGarrod fc Herbstl (120061 ) for hot cores and corinos, in which a gas-grain 
chemical model is utilized and both surface and gas-phase reactions are considered during the 
cold and heat-up phases. If L1527 is a representative example of a simple, isolated, collapsing 
core, perhaps more low-mass protostellar sources should exhibit a similar lukewarm envelope 
and the WCCC chemical signature. 
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Model 



We utilize the OSU gas-grain code, first described by lHasegawa. Herbst. fc Leung! (119921 ) . 
with some more recent modifications, the most important of which is the capability to in- 
crease the gas and dust temperatures with time ( Garrod fc Herbstlbood ). The model includes 
the granular parameters listed in Table [TJ We include the Rice-Rams perger-Kessel (RRK ) 
approach for reactive desorption from the solid phase, first suggested by lGarrod et al.l (120061 ). 
Reactive desorption, or non-thermal desorption of the products of exothermic surface reac- 
tions, is regulated by the parameter cxrrk, approximately equa l to the efficiency of the pro 
cess. For the models reported here cxrrk = 0.01, as chosen by lGarrod. Wakelam. fc Herbst 
( 120071 ). The results of corresponding models with orrk = are not significantly different 
from those presented, and are thus omitted for simplicity. In addition to the previously 
published modifications of the gas-grain network, we a dded the photodesorption of CO ice, 
based on the laboratory measurements of lOberg et al.l (120071 ). The desorption can be caused 
by both the ambient UV field and by secondary photons generated by cosmic rays. Because 
L1527 is a dense region, and not UV-dominated, we assume a typical interstellar radiation 
field and a value A v = 10 for all models presented. Consequently, the inclusion of CO 
photodesorption does not introduce a significant departure from results in which it is not 
included. The network includes 6312 reactions involving a total of 655 gaseous and surface 
species. Negative molecular ions have not yet been included in the model; an explanation of 
their abundances in L1527 is being undertaken using a purely gas-phase model. 

Our overall method is to follow the chemical evolution of a homogeneous (one-point) 
parcel of material as it progresses through the cold pre-stellar core phase and then heats 
up. It is assumed that the dust and gas temperatures are in equilibrium with one another. 



This a pproach differs from the hydrodynamic multi-point collapse model of lAikawa et al. 



(120081 ) in that we make simplified estimates of the time dependence and final value of the 
temperature and focus on the details of chemical evolution of groups of molecules during 
both the cold and warming-up phases. The advantage of our approach is that we are not 
wedded to a particular dynamical time dependence, the disadvantage is that it is unclear 
whether the physical changes that occur are totally realistic. 

We start the calculations by simulating the gas-grain chemistry that occurs in the central 
region of a cold pre-stellar core with a constant density of tir = 10 6 cm -3 and an initial 
temperature of T = 10 K. We start with gaseous atoms and ions except for molecular 
hydrogen, with abundances followi ng: low-metal elemental abundances as listed in Table [2] 
(IGarrod. Wakelam. fc Herbstl 120071 ). In addition, we consider the effect of an increased 
sulfur abundance of 1 x 10~ 5 , initially in the form of S + , and we discuss the results of this 
variation below. Figure [T] shows the fractional abundances of the major ice species H 2 0(s), 
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CO(s), and CH4(s) as functions of time under these conditions. The surface abundance 
of the predominant component, H 2 0(s), reaches a steady-state fractional abundance of ~ 
10~ 4 after a little more than 10 4 yr, while both CH 4 (s) and CO(s) reach a quasi-steady 
state at this time. After 10 6 yr, CO(s) decreases dramatically because it is hydrogenated, 
while CH4(s) actually increases in abundance. The initial cold phase in our warm-up models 
is allowed to continue for 10 5 yr, corresponding to the maximum CO(s) abundance, but 
we also run models to 10 4 and 10 6 yr for comparison. The abundance of CH 4 (s) in this 
range with respect to H 2 0(s) is calculated to be « 3.5%, which is in agreement with recent 
Spitzer Space Teles cope observ a tions according to SSH Y, and with previo us ISO results 



of 1-4% reported by iGibb et all (12004 ) . More recently, lOberg et all (120081 ) reported that 
the CH4(s) abundance may vary by 2-13% among a variety of sources, but that it is more 
narrowly defined to be 4.7 ± 1.6% if A^h 2 o(s) > 2 x 10 18 cm" 2 . Our calculation is in strong 
agreement with this new result. This extent of variation does not strongly affect the resulting 
abundances of the carbon-chain molecules to the best of our estimation. But, we cannot 
easily vary the amount of solid methane in our gas-grain model, since it is determined by the 
chemical history of the cold stage. In a gas-phase model in preparation (Harada & Herbst), 
we are looking into the dependence of assorted abundances on the gas-phase abundance of 
methane immediately after evaporation. 

After t he initial cold phase, the gas and dust temperatures are increased using the 
expression (IGarrod fc Herbstll2006l ) 



T = T + (T-a 



To) 



Aty 



where T is the initial temperature of 10 K, T max is the maximum temperature, is the 
heating timescale, and n is the order of the heating. If L1527 is in a transitional evolu- 
tionary stage between cold and hot phases, the temperature may be increasing above 30 
K. To account for this, the simulated temperature increase is carried to three maximum 
temperature values of T max = 30 K, 100 K, or 200 K in separate model runs. After reaching 
the asymptotic maximum, the model holds the temperature at that constant value while 
continuing to trace the chemical evolution. The adopted heating timescale is 2 x 10 5 yr with 
an order of n = 2 for each T ma ^, ba s ed on the estimated warm-up timescale for intermediate 



star formation. IGarrod fc Herbstl (120061 ) primarily utilized this timescale to simulate the 
chemistry during the formation of hot cores. 
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3. Results 

Our results are obtained as fractional abundances with respect to nn, Xi = rii/n^. 
The SSHY observations were reported as column densities, so comparison with observa- 
tion requir es us to convert these using the total h ydrogen column density iV H = 6 x 
10 22 cm" 2 ( jJOrgensen. Schoier. fc van Dishoeckl liool ). The converted observational values 



are listed in Table [3] along with the observed fractional abundances from the earlier study 



of iJorgensen. Schoier. &: van Dishoeckl (120041 ) and presented where appropriate in figures 



discussed below. The results of our simulations can be plotted as temporal evolutions of 
molecular abundances, both in the gas and on granular surfaces, both during the cold initial 
period and the subsequent warm-up, if it occurs. Two representative examples for unsatu- 
rated carbon chain gaseous species are shown in Figure [2J for C4H, and Figure [31 for HC5N. 
In these figures, we display the results for constant temperature, T = 10 K, and for heating 
to T max = 30 and 200 K. Models with similar heating but without surface reactions are 
also shown, but they do not show any significant distinction among temperature variations. 
These example plots also include the SSHY observations as solid horizontal lines and the 
temperature profiles of the warm-up models. The heating onset time was chosen to be 10 5 yr, 
as discussed above. Similar procedures that varied the heating onset time between 10 4 and 
10 6 yr showed no substantial difference in the results as a function of time after the start of 
heat-up. 

The dominant features for C4H and HC5N in models in which heating begins to occur 
after 10 5 yr include (i) a "cold" peak at 10 4-5 yr, which is essentially the "early-time" peak of 
cold gas-phase models, followed by a decrease due both to gas-phase reactions and accretion 
onto dust particles, (ii) an increase in abundance that occurs when volatile molecules such 
methane and carbon monoxide begin to sublimate from the grain surfaces as the temperatures 
nears 30 K, and (iii) a second, or warm, peak, possibly with fine structure, which occurs in 
the time frame 10 5-6 yr, and (iv) a final decrease in abundance. The second peak occurs 
at or near the period when the maximum temperature in the model (30 K or 200 K) is 
reached. The warm peak abundance for C4H and HC5N is larger for the higher temperature 
model; in addition, the abundances when the higher temperature model is temporarily at 
30 K are very similar to the abundances achieved when the 30 K warm-up model reaches its 
asymptotic temperature. In the case where T max = 200 K, a satellite warm-up peak occurs 
near T = 30 K and has abundance similar to that of the single T max = 30 K warm-up peak. 
The second, larger peak occurs when the temperature first reaches 200 K and the majority 
of the ice mantle species are driven from the grains. 

In the cold model, the abundances of C4H and HC5N are slightly regenerated in a plateau 
region beginning at t ~ 6 x 10 5 yr until t ~ 4 x 10 6 yr as a result of reactive desorption. 
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Actually, the regenerated abundances of most molecules do not exceed the initial cold peak 
with the exceptions of Ci , C 5 N, CS, H 2 CS, and H 2 S. The species CN and HCO + do not 
decrease from the initial cold peak but remain at a plateau value for nearly 8 x 10 6 yr. Finally, 
as exemplified by C4H and HC5N, the warm up in the absence of surface chemistry shows 
no warm peak, demonstrating that surface chemistry prior to warm-up is indeed necessary 
for the production of gas-phase molecules in the vicinity of 30 K. The detailed mechanism 
for the production of the gas-phase molecules during the warm-up is discussed below. 

In Table |3] we list peak abundances for a large number of species for four models in 
addition to observed fractional abundances for L1527 if available. There are three columns 
that specify a maximum asymptotic temperature of 30 K, 100 K, and 200 K. Here the peak 
values refer always to that achieved during warm up after t = 10 5 yr and not in the cold 
phase. In models with an asymptotic temperature that exceeds 30 K, the peak values can 
be greater than values in the vicinity of 30 K, as is the case for most unsaturated carbon 
chain species (e.g., C 4 H and HC 5 N) as well as all standard hot-core- type molecules. The 
potentially controversial high-temperature results for the unsaturated species are discussed 
in Section 3.3. In the column in which results from a model with a constant temperature of 
10 K are listed, the peak abundances are the cold, "early-time" values prior to t = 10 5 yr. 
One can see that for unsaturated classes of organic molecules such as the radicals C n H, and 
the cyanopolyynes HC2 n CN, the peak abundances during the warming up periods are larger 
than for the cold period. On the other hand, for other classes of species (e.g. molecular ions) 
the issue is less clear. 



It is desirable to understand the mechanism(s) responsible for the increase in abundance 
of gas-phase carbon-chain species during the warm-up period. As gas-phase species freeze 
out onto grains in the cold era, they are regenerated when highly volatile precursor species 
sublimate as the temperature approaches ~ 30 K and subsequently react in the gas phase. 
Methane plays a key role in the formation of unsaturated hydrocarbons and so its evolution 
is shown in Figure H] for models with no warm-up and warm-up to 30 K and 200 K. During 
the cold era, a large abundance of methane ice builds up in all models, formed mainly by 
sequential hydrogenation of carbon on surfaces, in which the concluding reaction is 



The desorption of methane that occurs as a result of heating to 30 K and beyond becomes 
apparent in the upper panels of Figure |U As this desorption occurs, gas-phase chemistry 



3.1. Chemical Interpretation 




(2) 
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begins to regenerate many unsaturated species. The corresponding increase in C4H abun- 
dance is evident in Figure [2] at t ~ 3 x 10 5 yr for T max = 30 K and t ~ 1.5 x 10 5 yr for 
T max = 200 K. Models without grain-surface reactions form much less methane (see the 
abundances of CH 4 (g.o.) in Figure H|), and thus are incapable of regenerating secondary 
products, as can be seen for C4H and HC5N in Figures [2] and [3] in the curves labeled "g" . 

The regeneration of unsaturated carbon-chain molecules does not occur as simple sub- 
limation from grain surfaces, but rather as the result of subsequent gas-phase chemistry, simi- 



l ar to that discussed in a dynamic model of the TMC-1 sub-cores by lMarkwick. Millar, fc Charnley 



( 120001 ). For the models with T max = 30 K the most fundamental steps following the subli- 
mation of methane are those leading to the formation of acetylene and related hydrocarbon 
ions: 

C+ + CH 4 — > C 2 H+ + H, (3) 
C + + CH 4 — ► C 2 H+ + H 2 , (4) 

followed by 

C 2 H+ + H 2 — > C 2 H+, (5) 
which lead to the dissociative electron recombination reactions 

C 2 H+ + e~ — ► C 2 H 2 + H (6) 

and 

C 2 H+ + e" — > C 2 H 2 + 2H. (7) 

Further reactions among the two-carbon ions and acetylene produce three- and four-carbon 
atom species. These products react in turn to form the longer chain C n H and C n H 2 species. 
A by-product of the ion molecule network, CN will also react with acetylene: 

CN + C 2 H 2 — > HC 3 N + H, (8) 

which is the primary formation path of the first member of the polyyne family. Further, 
more complex reaction pathways form longer C n N and HC 2n CN species. 

In models where the temperature increases to 100-200 K, t he unsaturated ch e mistr y 



discussed here is combined with the hot-core chemistry discussed in lGarrod fc Herbstl (120061 ). 
in which a radical-based surface chemistry followed by the desorption of less vaporizable 
molecules helps to produce the oxygen-rich organic molecules detected in hot cores. The 
observed and peak calculated abundances of some typical hot core molecules such as CH 3 OH, 
H 2 CO, and CH3CN are included in Table HI where it can be seen that the lukewarm heating 
is insufficient to reproduce the observations of typical hot cores. 
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Our warm-up models treat the desorption of surface species as if they were pure ices. 
In reality, mixtures of ices act in a more complex manner. The surface experiments of 



Collings et al.l (120041 ) address the trapping of more volatile species where they are either 



deposited as a surface layer on top of amorphous water ice, or co-deposited with water. 
The results indicate that a significant fraction of the most volatile molecules, CO, N2, and 
O2, can be partially retained in the icy mantle to temperatures much greater than their 
ty pical sublimation temperatures. Such an effect has been considered in the hot-core models 
of IViti et al.l (120041 ) . Our gas-grain code currently does not model surfaces in sufficient 
detail to account for this retentive effect of water; however, comparison with the other 
molecules would indicate that not more than half of the available CH 4 (s) would likely be 
retained. The fraction released to the gas phase under these conditions still represents a 
tremendous enhancement over the amount of CH 4 (g) formed by gas phase processes only. 
Moreover, a gas-phase model in preparation by Harada & Herbst shows that the production 
of unsaturated species from methane is not strongly affected by such a diminution in the 
amount of gaseous methane vaporized in the vicinity of 30 K. 



3.2. Comparison with observations 



Figures [2] and [3] have an important contrast that helps to clarify the major question 
of this paper. The cold "early-time" peak of C4H is less than the observation by more 
than an order of magnitude, and only the warm-up models produce enough C4H to agree 
with the observation. However, the observed abundance of HC5N can be fit to within an 
order of magnitude with warm-up models or with the cold early-time peak. We therefore 
compare the results of all the observed species to determine if the composition of L1527 
requires heating to match the observed composition, or if the results are ambivalent. The 
reaction network does not differentiate between the normal and carbene isomers of the C n H2 
molecules, which have been detected in L1527. We assume that the observable vinylidene 
isomer H 2 C 3 represents 2% of the total abundance of C3H?, as calculated for a dark cloud 
using the gas-phase network described in lPark. Wakelam. fc Herbstl (120061 ). We extend this 
estimate to the entire group of C n H 2 species in order to compare the calculations with the 
abundances of isomers that were actually observed. 

To formulate a quantitative answ er to the cold-warm up dichotomy, we revisit the "mean 
confidence level" method utilized by iGarrod. Wakelam. fc Herbstl ( 120071 ). The confidence 
level, Ki, for agreement between the computed abundance, Xi, and the observation, X b s ,i, 
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for each relevant molecule is computed as 

where the agreement improves as Ki approaches 1. As before, we define o = 1 so that 
Ki = 0.317 indicates a one order of magnitude departure from the observed value. In addition 
to the SSHY observations, we include the JSV04 observations to analyze the composition of 
the cloud. The average value k of the confidence level for each set of observations and for 
all observations combined is shown in Figure [5] as a function of time. 

Within this figure, several meaningful peaks arise. Prior to the heating, the mean 
confidence level peaks at early-time for all three data sets in the region of the cold peak. For 
the heated models, a peak appears at t ~ 3 x 10 5 yr for T max = 30 K and at t ~ 1.5 x 10 5 yr 
for T max = 200 K. These peaks indicate that the cold models show an overall agreement 
within one order of magnitude of the observation, and that the heated models also agree 
with observation to a similar extent when T pa 30 K. There is also a secondary peak for 
the 200 K warm-up models at higher temperatures, which indicates that abundances have 
increased to overly large values compared with observation as the temperature exceeds 30 
K, and then decrease once again. 

Now let us compare the data sets more carefully. The bottom panel in Figure [5] displays 
the SSHY observations of unsaturated carbon-chain and cyanopolyyne molecules. Here it 
can be seen that the confidence level peaks associated with heating exceed that in the cold 
period. Specifically, the confidence level associated with heating is elevated to 0.628 at 
t = 2.7 x 10 5 yr and T = 25 K for the T max = 30 K model. The peak confidence level for 
the T max = 200 K model is 0.663 at t = 1.6 x 10 5 yr and T = 27 K, indicating average 
agreement for heated models to within a factor of 3. The confidence level of the cold peak 
is 0.380 at t — 3.5 x 10 4 yr, indicating agreement to within a factor of 7.5. In this set 
of observations, heating gives an overall better fit to the observed data. The middle panel 
displays the JSV04 observation set. This plot shows an ambivalence between the cold peak 
and the heated models, as the maximum confidence levels are 0.631 for the cold peak at 
t = 6.5 x 10 4 yr, 0.648 for T max = 30 K at t = 2.9 x 10 5 yr and T = 29 K, and 0.606 for 
T m ax = 200 K and T = 27 K at t = 1.6 x 10 5 yr, indicating agreement to within a factor of 3 
for all three possibilities. The two sets of observations are combined in the top panel, where 
the warmed models show slightly higher confidence levels of 0.624 and 0.639 as compared 
with 0.476 (factor of 5) for cold processes. The optimal agreement peaks for the combined 
data set occur at the same times as for the SSHY observations. 

This analysis marginally favors the explanation that the observed L1527 composition 
results at least partially from heating, although it should be noted that the very high confi- 
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dence level for the warm-up peak lasts only a brief period of time (~ 10 4 yr), especially for 
models with higher maximum temperatures, where the post-peak drop is precipitous. The 
significance of dividing this confidence analysis by observation set is that the SSHY results 
included a profile mapping procedure to demonstrate that the velocity of C4H observations is 
associated with the central protostar and spread throughout the gaseous envelope of L1527, 
rather than associated with the outflows. The JSV04 observations do not distinguish sys- 
tematic velocity components. Future study of this region with detailed interferometry would 
certainly add more spatial abundance information and likely better resolve the dichotomy 
between cold and warm explanations for the composition. 

Some individual molecular abundances are summarized in Figures |6] and [7J at the opti- 
mum times determined from the above analysis. Figure El compares the optimal warm-up 
results for T max = 30 K with the optimal cold results, whereas Figure [7] is a similar com- 
parison for the optimal T max = 200 K results. The optimal cold results pertain to a cold 
pre-stellar source with density 10 6 cm -3 and should not be confused with dark cloud abun- 
dances such as those of TMC-1. The bottom panels of these figures display the results for 
the SSHY species that do not require any adjustment, with the observed values denoted by 
asterisks where appropriate. The center panels include the C n H 2 species detected by SSHY 
and some potentially detectable varieties of carbon-chain molecules. The top panels display 
the results relating to observations of JSV04. In general, these plots also indicate on a visual 
basis that the SSHY molecules are well reproduced by warm-up models. Optimal abun- 
dances are also listed in Table which includes the observed molecules as well as a larger 
number of predictions. In this table, italic type indicates that the computed abundance 
exceeds the observed value by one order of magnitude or more. Boldface type indicates that 
the computed abundance is too small by at least one order of magnitude. 

From Table \5\ we can see the individual levels of agreement for observed species at 
optimum times. At the optimum time for agreement of warm-up models, a total of 17 of 22 
of the observed molecules are within one order of magnitude of the observation, including 
11 of the 13 SSHY species and 6 of the 9 from JSV04. Variations in heating do change the 
optimum time but not the number of molecules in agreement. The molecules which do not 
agree with observation at the optimum time are C5H, CgH, CS, CN, and N 2 H + , where the 
models over-produce the first four, but do not make enough N 2 H + . 

At the optimum time of agreement for the cold-phase model, t = 3.5 x 10 4 yr, only 5 of 
the 13 SSHY molecules agree with observation to within an order of magnitude, while 8 of 
9 JSV04 molecules agree, for a total of 13 out of the 22. One reason for this comparatively 
worse agreement can be seen in Figures [2] and [31 where the cold peaks arise at different times. 
Overall, cold processes in this model tend to under-produce the carbon chain molecules 
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observed by SSHY, and cannot produce C4H or HCgN to within an order of magnitude 
of the observation at any time. The superiority of the warm-up model for carbon-chain 
species is seen dramatically in Figures M and [7J As the carbon chain molecules get longer, 
the abundances relative to the cold values are drastically enhanced. This effect occurs to 
such an extent that longer chain molecules of the classes C n H, C n H2, and HC2 n CN have 
only slighter lower abundances than the shorter chains in that group, at least on the scale 
depicted here. The trend presents a good criterion for prediction of further observations of 
lukewarm corinos. As shown in Figure 6 of SSHY, the abundances of carbon-chain species 
observed in L1527 tend to be slightly lower than those observed in the cold core TMC-1. 



These latter o bservations are fit by cold models in which car 



must be used ( ISmith. Herbst. fc Chang 



thermal mantle evaporation occurs (IMarkwick. Millar. Sz Charnleyll2000l ) 



Don-rich elemental abundances 



200 i 



Wakelam et al 



2006), or m odels in which non- 



Among the molecules detected in L1527, two stand out as almost uniquely abundant: 
C2H and C3H2. SSHY demonstrated that the abundance of C2H is remarkably high in 
compari son with seyeral o ther star-forming and starless cores on the basis of a previous 
survey (jSuzuki et al.lll992l ). The calculated abundance of C2H for T max = 30 K remains 
elevated at this temperature at nearly one order of magnitude greater than the cold peak 
value for several million years. However, the abundance of the T max = 200 K model peaks 
and then rapidly declines to the level of the cold value as T approaches 100 K, prior to 
t = 2 x 10 5 yr. This unusual result explains that L1527 at T w 30 K can reasonably have 
a C2H abundance that is nearly an order of magnitude greater than both colder starless 
cores and hotter sta r forming regions. This finding is in agreement with the C2H image of 
Beuther et al.l (120081 ) for IRAS 18089-1732, which shows a hole toward the hot core. 



Also, iBenson. Caselli. fc Myersl (119981 ) indicated that the abundance of C3H2 was re- 
markably greater toward L1527 and the TMC-1 sub-cores than toward any other cold sources 
included in their survey. Our calculated lukewarm abundance of C3H2 increases by an order of 
magnitude over the cold abundance and remains elevated for a large period of time. The high 
abundances of both C2H and C3H2 lend support to the explanation that partial mantle re- 
moval in L1527 results in system atic enhancements in carbon-chain abundances, although th e 
case is less obvious for TMC-1 ( IWakelam et al.ll2006l ; IMarkwick. Millar, fc Charnleyll2000l ) 



3.3. Observational Predictions 

The peak abundance values for a variety of carbon-chain and other molecules for T max = 
30 K are listed in Table [3] while the abundances at the optimal times for L1527 are listed 
in Table for all warm-up models. This information provides observational predictions for 
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L1527, as well as any other potential lukewarm corino sources. 



While our model calculations indicate that warming can successfully regenerate carbon- 
chain species, they also introduce an unexpected result for hot-core and hot-corino chemistry. 
Heating to T max = 200 K increases the peak abundances of these molecules by at least an 
order of magnitude above the T max = 30 K peak (see Table [3]). Moreover, except for the 
species C2H (see above), these high abundances do not decline rapidly. Previously, carbon- 
chain spec ies have not been reported in hot core /corino observation s with the exception 
of HC 3 N JWvrowski et al.l Il999l : bibb et al.l I2OO0I ; felake et all Il987h . to the best of our 
knowledge. It is not clear if unsaturated chain hydrocarbons are simply not present in the 
observed hot cores, or if they are not detected in the complex spectra of these regions, 
which are typically dominated by so-called "weeds". However, in a recent survey of the 
hot corino source IRAS 16293+2422, a number of small carbon-chain species have been 
detected and are likely to exist in the warm region of the source (E. Caux, C. Ceccarelli, 
private communications). This survey lends a degree of support toward our model result of 
enhanced carbon-chain abundance in hot regions. However, it is also possible that the high 
abundances we predict in hot cores may be the result of a lack of destruction mechanisms 
that become important at higher temperatures. 



3.4. Increased Sulfur 

The gas-phase elemental abundance of sulfur in cold dense sources is highly uncertain 
because there is no evidence for depletion in diffuse clouds, yet use of a high value for cold 



dense sources tends to worsen agreement with observation (IQuan et al.ll2008l ). If the sulfur 
abundance in our model is raised to 1 x 10" 5 (see Table certain results are affected much 
more strongly than others. For the cold phase, the major changes are large diminutions in 
the peak abundances of the C n H species and an increase of at least an order of magnitude 
for the sulfur-bearing species such as CCS. For the warm-up phase of the T max = 30 K 
model to lukewarm corino conditions, the results change little, with peak abundances for all 
species listed in our tables altered at most by a factor of three. For the T max = 100 K and 
7m ax = 200 K models, the peak abundances tend to be reduced for the carbon-chain species, 
and their post-peak declines tend to occur more quickly, at a time scale of at most 10 6 yr, 
reducing the possible discrepancy with observation for carbon-chain species in hot cores of 
this age. 

With the mean confidence level approach, we see that the optimum n for the observed 
species in the cold era decreases to 0.350, barely corresponding to order-of-magnitude agree- 
ment. The optimum k for the warm-up models does not change significantly, on the other 
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hand, indicating a stronger preference for lukewarm processes than in the low-sulfur case. 



4. Conclusions 

Using our gas-grain chemical model, we have attempted to distinguish the chemical 
composition of L1527 from more familiar evolutionary phases. More specifically, this com- 
position poses the dilemma that the observed gas-phase abundances might have required 
heating, might be explained by cold processes alone, or might require a combination of the 
two depending upon exactly where the molecules are located within the source. We have 
approached this question by computing a series of one-point model simulations of a dense 
core that include various degrees of heating after a preliminary cold phase lasting 10 5 yr. 
A mean confidence level analysis indicates that the composition of the cloud is somewhat 
more likely to be a result of heating and regeneration than of cold processes. Dividing the 
observed molecules into two classes - the carbon-chain species detected by SSHY and the 
species detected by JSV04 - we find that the preference for a warm-up model is stronger 
for the former species although the level of high agreement during the warm-up period is 
brief. This preference is strengthened by the use of a high value for the sulfur elemental 
abundance. 

The best agreement between model and observation occurs as the most volatile molecules 
are being desorbed from the grain surfaces. This agreement occurs at t — 2.7 x 10 5 yr and 
T = 25 K for the T max = 30 K model and t = 1.6 x 10 5 yr and T = 27 K for the T max = 200 K 
model. Of the volatile species, the most important chemically is CH4 because it acts as a 
precursor for a gas-phase chemistry leading to carbon-chain species, as discussed by SSHY . 



The distinctly large abundances of C 2 H (SSHY) and C 3 H 2 ( iBenson. Caselli. fc Myers! Il998l ) 
as compared with various other sources provide some additional support for this lukewarm 
explanation of the L1527 composition. Further, the gradual drop off in abundance of carbon- 
chain molecules in L1527 with increasing size is fit better with our warm-up than our cold 
model. It should be noted, however, that our model here uses oxygen-rich elemental abun- 
dances, and that such models cannot reproduce the abundances of complex molecules seen 
in TMC-1 , so may underestima te the extent of molecular complexity under cold conditions 



in L1527 (IWakelam et al.ll2006l ). Although the envelope of L1527 at 30 K is warmer than a 
cold core at 10 K, it is not nearly as hot as a hot core or hot corino. The chemistry is clearly 
distinguishable from the more evolved phase, because the temperature is not high enough 
to release some products of grain-surface chemistry such as CH3OH. If L1527 is typical of 
a borderline Class / I low mass protostar, then the "lukewarm" corino phase could be 
fairly common. Observations of L483 show similar evolutionary characteristics including an 
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enhanced dust temperature of Td = 40 K (ITafalla et al.ll2000i ); thus, this source may be a 
good candidate for a molecular survey of warm carbon chain chemistry. 

One complication that arises is that the models in which the temperature rises to 100 
K or 200 K produce large abundances of most gaseous carbon-chain species in addition to 
the standard hot-core species and these abundances linger on timescales of 2 — 5 x 10 6 yr 
after enhancement, unless a high elemental abundance for sulfur is utilized. The most salient 
exception to this rule is C2H, for which we predict a low abundance in hot cores and corinos. 
The prob lem of high ca r bon-c hain abundances in hot cores also occurs in the hydrodynamic 
model of lAikawa et al.l (120081 ). in which large abundances of these species are maintained 
as the gas collapses to the central star and warms up. The peak results from Table [3] show 
that these abundances are somewhat lower than abundances for previously reported hot core 
molecules (see Table HI), but may be large enough to be detectable unless their widely-spaced 
spectral lines are hidden by those of weeds. Unfortunatel y, only detections of HC3N toward 



hot cores have previou sly been reported in the literature (IWyrowski et al.lll999l ; iBlake et al. 



19871 ; iGibb et al.ll2000l ) . Some of these detections are sufficiently strong that one wonders why 
more complex cyanopolyynes have not yet been detected. Low upper limits towards Orion 
have been reported to us by several investigators but this source is complex and may not be 
representative. A report of a recent survey in progress (E. Caux, private communication) 
indicates, on the other hand, that some smaller carbon-chain species are indeed observable 
in the hot corino IRAS 16293+2422, most likely associated with the warm region of the 
source. Further observations of carbon-chain molecules, particularly with lengths of four 
or more carbon atoms, toward additional hot core or corino sources might strengthen our 
characterization of LI 52 7. 

E.H. would like to thank the National Science Foundation for supporting his research 
program in astrochemistry and NASA for support to study interstellar ice formation. R.T.G. 
thanks the Alexander von Humboldt Foundation for a research fellowship. 
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Fig. 1. — Surface fractional abundances of major mantle species with respect to total hydro- 
gen at rin = 10 6 cm -3 and T = 10 K. 
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Fig. 2. — Temporal evolution of the C4H abundance. All models undergo a cold pre-stellar 
core period of 10 5 yr. The results are computed for three temperature variations: constant 
T = 10 K, labeled (a); warm-up to T max = 30 K, labeled (b); warm-up to T max = 200 K, 
labeled (c). Models without grain surface chemistry are included (g). These models overlap 
for all three heating variations. The temperature profile for the two warm-up models is also 
plotted vs time. The observed abundance is indicated with a horizontal line. 
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Fig. 3. — Temporal evolution of the HC5N abundance. The labels are the same as those 
used for Fig. [2j 
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Fig. 4. — Abundance evolution of CH4 with heating. Solid lines represent CH 4 (s) and dashed 
lines represent CH 4 (g) in gas-grain models. For comparison, dotted lines labeled CH 4 (g.o.) 
represent CH 4 (g) formed by gas-phase processes only. Panel (a) contains results for constant 
T = 10 K, panel (b) for warm-up to T max = 30 K, including the temperature profile, and 
panel (c) for warm-up to T max = 200 K, also with a temperature profile. 
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Fig. 5. — Mean confidence level k for observe d molecules i n L152 7 as a function of time. 
Results for 13 carbo n-chain species observed by ISakai et all (I2008af) in panel (a), for 9 other 



species observed by IJorgensen. Schoier. fc van Dishoeckl (120041 ) in panel (b), and for all 22 
observed species in panel (c). The solid line represents T max = 30 K models, the dashed line 
represents T max = 200 K models, and the dotted line represents cold models. The horizontal 
line represents one order of magnitude average difference between observed and calculated 
abundances. 
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Fig. 6. — Calculated abundances of various carbon-chain and other molecules at the time 
of optimum agreement between models and observations for warm-up to T max = 30 K. 
Abundance results for T max = 30 K (□) at t = 2.7 x 10 5 yr are compared with optimal cold 
values for constant T = 10 K (■) at t — 3.5 x 10 4 yr. Observations (*) are denoted where 
appropriate. 
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Table 1. Granular Properties 



Parameter 




value 




Sticking coefficient 


S 


0.5 




Radius 


a 


0.1 


fim 


Density 


P 


3 


g cm -3 


Mass Fraction 


x d 


0.01 




Reactive Desorption Efficiency 


«RRK 


0.01 
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Table 2. Initial Abundances 



Species i rii/nn t 



H 2 


0.5 


He 


0.14 


C+ 


7.3(-5) 


N 


2.14(-5) 





1.76(-4) 


s+ 


8.0(-8) a 


Na+ 


2.0(-9) 


Mg+ 


7.0(-9) 


Si+ 


8.0(-9) 


P+ 


3.0(-9) 


C1+ 


4.0(-9) 


Fe+ 


3.0(-9) 



to(6) = a x 10 b 

a We also consider the 
effect of S+=1.0(-5) 



Table 3. LI 527 & Peak Calculated Fractional Abundances of Assorted Species 



k.j l_. L. 11... o 


Obs. 


T = 10 K 


T — 3D K 

± max — iv 


T — 1 no K 

- 1 max — iuvj 


T — 200 K 

± max — £ uu xv 


Ohq Rpf 


C 2 H 


8.3(-09) 


2.0(-09) 


8.9(-09) 


1.0(-08) 


l.l(-08) 


a 


C 3 H 




5.0(-10) 


1.9(-09) 


1.2(-08) 


1.5(-08) 




C 4 H 


3.3(-09) 


1.8(-10) 


3.0(-09) 


3.3(-08) 


4.1(-08) 


a 


C 5 H 


L6(-ll) 


1.2(-11) 


8.7(-10) 


8.1(-09) 


2.9(-08) 


a 


C 6 H 


l.O(-ll) 


1.5(-11) 


6.0(-10) 


5.8(-09) 


3.5(-08) 


a 


C 7 H 




9.3(-13) 


4.4(-10) 


2.7(-09) 


3.6(-08) 




C 8 H 




1.2(-12) 


7.2(-10) 


3.1(-09) 


3.1(-08) 




C 9 H 




4.8(-14) 


5.1(-10) 


2.4(-09) 


3.4(-08) 




HC 5 N 


9.7(-ll) 


9.1(-11) 


5.4(-10) 


5.1(-08) 


2.0(-07) 


a 


HC 7 N 


2.7(-ll) 


4.3(-12) 


1.9(-10) 


6.3(-09) 


4.1(-08) 


a. 


HCgN 


2.5(-12) 


1.5(-13) 


9.9(-ll) 


1.7(-09) 


1.7(-08) 


a 


C 2 H 2 




6.5(-09) 


4.2(-08) 


1.6(-07) 


2.3(-07) 




C3H2 


l.O(-ll) 


9.5(-13) 


1.8(-11) 


4.8(-10) 


4.7(-10) 


a,b 


C4H2 


2.2(-ll) 


2.9(-12) 


1.7(-10) 


1.4(-09) 


2.2(-09) 


a , b 


C5H2 




1.4(-13) 


4.8(-ll) 


1.3(-09) 


1.4(-09) 


b 


C 6 H 2 


2.5(-12) 


1.2(-13) 


6.2(-ll) 


2.0(-10) 


1.0(-09) 


a,b 


C7H2 




9.9(-15) 


2.1(-11) 


1.4(-10) 


6.0(-10) 


b 


C 8 H 2 




8.7(-15) 


3.0(-ll) 


7.1(-11) 


4.0(-10) 


b 


C9H2 




5.9(-16) 


2.0(-ll) 


9.7(-ll) 


l.l(-09) 


b 


CH3CCH 


7.8(-10) 


2.4(-10) 


4.0(-09) 


9.6(-07) 


9.3(-07) 


a 


CH3C4H 




2.0(-12) 


2.3(-09) 


7.2(-09) 


5.9(-08) 




CH 3 C 6 H 




3.1(-13) 


5.8(-10) 


2.1(-09) 


7.6(-09) 




CCS 


8.5(-ll) 


2.0(-10) 


1.3(-10) 


1.2 (-09) 


1.1 (-09) 


a 


HCO^ 


1.0(-12) 


1.6(-12) 


7.8(-13) 


2.4(-ll) 


3.8(-12) 


c 


HCN 


1.2(-09) 


2.1(-07) 


6.9(-09) 


1.7(-06) 


1.8(-06) 


d 


HNC 


3.2(-10) 


1.4(-07) 


5.6(-09) 


7.7(-07) 


1.2 (-06) 


d 


HC3N 


8.9(-10) 


3.1(-09) 


8.8(-10) 


1.2(-07) 


1.4(-07) 


d 


CO 


3.9(-05) 


3.3(-05) 


5.2(-05) 


7.3(-05) 


7.3(-05) 


d 


CS 


3.3(-10) 


6.4(-08) 


8.0(-09) 


5.7(-08) 


5.0(-08) 


d 


so 


1.4(-10) 


1.1 (-09) 


1.5(-10) 


5.6(-08) 


5.0(-08) 


d 


HCO+ 


6.0(-10) 


4.4(-10) 


1.8(-09) 


2.2(-09) 


1.8(-09) 


d 


N 2 H+ 


2.5(-10) 


6.3(-ll) 


6.1(-10) 


6.2(-ll) 


8.1(-11) 


d 


CN 


8.0(-ll) 


3.2(-09) 


4.0(-09) 


3.0(-07) 


3.0(-07) 


d 


c 2 




3.4(-09) 


2.0(-09) 


4.3(-08) 


3.6(-08) 




C 3 




2.3(-06) 


3.9(-08) 


2.7(-07) 


2.0(-07) 




c 4 




8.2(-09) 


1.3(-09) 


1.2 (-08) 


1.8(-08) 




Cg 




1.4(-08) 


1.5(-08) 


6.2(-08) 


2.4(-07) 




C e 




6.0(-ll) 


3.3(-09) 


5.3(-09) 


5.8(-08) 




C 7 




6.4(-10) 


7.8(-09) 


6.9(-09) 


1.7(-07) 




C 8 




5.8(-12) 


2.9(-09) 


2.3(-09) 


5.1(-08) 




Cg 




4.4(-ll) 


1.4(-08) 


7.4(-09) 


2.7(-07) 




C10 




8.2(-14) 


9.4(-10) 


1.4(-09) 


1.5(-08) 




C 2 N 




1.2(-10) 


8.9(-10) 


3.3(-08) 


4.5(-08) 




C 3 N 




2.2(-ll) 


3.0(-10) 


3.4(-08) 


3.3(-08) 




C 4 N 




1.6(-11) 


9.2(-10) 


3.2(-08) 


7. 2 (-08) 




C 5 N 




6.8(-13) 


5.6(-10) 


2.7(-08) 


6.6(-08) 
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Tablc 3 — Continued 



Species Obs. 


T = 


10 K 


Tinax — 


= 30 K 


Tmax — 


= 100 K 


Tmax — 


= 200 K Obs. Ref. 


C 7 N 


4.5 


-14) 


2.2( 


-10) 


3.0 


'-09) 


2.0( 


-08) 


C 9 N 


2.8 


-15) 


1.1( 


-10) 


4.8 


;-io) 


7.3( 


-09) 


CH3C3N 


7.5 


-12) 


7.3( 


-13) 


1.9 


;-n) 


1.8( 


-10) 


CH3C5N 


2.1 


-13) 


8.8( 


-14) 


2.1 


:-n) 


2.7( 


-10) 


CH3C7N 


9.5 


-15) 


3.2( 


-14) 


2.7 


;-i2) 


3.2( 


-11) 


N 2 


1.9 


-06) 


5.6( 


-06) 


1.1 


'-05) 


1.1( 


-05) 


2 


3.7 


-07) 


3.7( 


-08) 


4.8 


-05) 


4.3( 


-05) 


C0 2 


5.4 


-07) 


8.6( 


-09) 


9.2 


;-06) 


1.8( 


-06) 


C3H3 


2.1 


-10) 


3.9( 


-09) 


7.5 


;-08) 


8.5( 


-08) 


C4H4 


3.7 


-11) 


6.1( 


-10) 


2.6 


[-07) 


2.3( 


-07) 




5.4 


-13) 


1.6( 


-09) 


6.6 


[-09) 


1.8( 


-08) 


H 2 CS 


1.9 


-10) 


5.7( 


-10) 


1.8 


;-08) 


4.6( 


-08) 


H 2 S 


1.5 


-09) 


4.7( 


-09) 


8.7 


1-09) 


7.9( 


-09) 


CH 2 CN 


1.9 


-10) 


1.0( 


-10) 


6.4 


[-08) 


8.4( 


-08) 


CH 2 CO 


6.4 


-09) 


1.0( 


-09) 


1.3 


-07) 


1.4( 


-07) 



a reported in lSakai etaH l|2008al ) 

b calculations adjusted for linear isomer as described in § 13.21 
c unpublished, recently observed by N. Sakai (private communication) 
d reported in |j0rgensen. Schoier. 8z van Dishoeckl ll2004h 
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Table 4. Observed & Peak Calculated Fractional Abundances of Assorted Hot 

Core/Corino Species 



Species 


Obs. 


T = 10 K 


Tmax — 


: 30 K 


Tmax — 


100 K 


T max = 200 K 


Source 


HCN 


l(-06) a 


2.1 (-07) 


6.9(- 


■09) 


1.7(- 


■06) 


1.8(-06) 


G327.3-0.6 


CH 3 CN 


7(-07) a 
[1.6±0.2](-09) 6 
[1.0±0.4](-09) c 


3.4(-09) 


1.5(- 


■11) 


7.4(- 


■08) 


8.3(-08) 


G327.3-0.6 
NGC 1333 IRAS 4A e 
IRAS 16293 2422 e 


CH3OH 


2(-05) a 
< 7(-09) b 
3(-07) c 


6.0(-10) 


6.9(- 


■11) 


2.5(- 


■07) 


1.2(-06) 


G327.3-06 
NGC 1333 IRAS 4A e 
IRAS 16293 2422 e 


CH3OCH3 


< 2.8(-08) b 
[2.4±3.7](-07) c 


7.7(-15) 


3.8(- 


■16) 


4.2(- 


■10) 


3.2(-09) 


NGC 1333 IRAS 4A e 
IRAS 16293 2422 e 


H 2 CO 


2(-08) d 
l(-07) c 


1.2(-09) 


4.3(- 


■09) 


3.8(- 


■08) 


3.7(-08) 


NGC 1333 IRAS 4A e 
IRAS 16293 2422 e 


HC 3 N 


3(-ll) a 
1.6(-09)f 


3.1(-09) 


8.8(- 


■10) 


1.2(- 


■07) 


1.4(-07) 


G327.3-0.6 
Orion KL 


HCOOH 


[4.6±7.9](-09) 6 
6.2(-08) c 


1.8(-10) 


6.0(- 


■12) 


5.1(- 


■09) 


1.6(-08) 


NGC 1333 IRAS 4A e 
IRAS 16293 2422 e 


HCOOCH3 


[3.4±1.7](-08) b W 
[3.6±0.7](-08) 6 (") 
[1.7±0.7](-07) c « 
[2.3±0.8](-07) c <") 


1.5(-15) 


3.1(- 


■15) 


8.3(- 


■11) 


4.7(-10) 


NGC 1333 IRAS 4A e 
NGC 1333 IRAS 4A e 
IRAS 16293 2422 e 
IRAS 16293 2422 e 


HNCO 


3(-09) a 


4.6(-09) 


4.5(- 


12) 


7.6(- 


■07) 


8.3(-07) 


G327.3-0.6 



HGibb et al] bOOOh 

iBottinelli et al.1 J2004h : (i) A substate of HCOOCH3 (ii) E substate of HCOOCH3 
fcazaux et all J2003h : (i) A substate of HCOOCH3 (ii) E substate of HCOOCH3 
^Maret et ail II2004TI 
e hot corino 
iBlake et all Jl987h 
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Table 5. Observed and Optimal Time Fractional Abundances 







max — iv 


max — xv 


T — 

J max — 


100 K 


max — ^yJKJ xv 


k. j yj l\_.o 


Obs. a 


X 


X 




opt 


X 


C 2 H 


8.3(-09) 


3.9(-10) 


8.9(-09) 


1.0( 


-08) 


1.1 (-08) 


C 3 H 




4.9(-10) 


1.0(-09) 


8.6( 


-10) 


8.4(-10) 


C 4 H 


3.3(-09) 


8.2(-ll) b 


1.7(-09) 


1.1' 


-09) 


9.9(-10) 


C 5 H 


1.6(-11) 


1.2(-11) 


4.6(-10) c 


3.3 


'-10) 


3.1 (-10) 


C 6 H 


l.O(-ll) 


1.3(-11) 


2.3(-10) 


1.3 


'-10) 


1.1(-10) 


C 7 H 




7.7(-13) 


1.3(-10) 


6.9 


-11) 


6.1(-11) 


C 8 H 




5.9(-13) 


l.l(-lO) 


4.5( 


-11) 


3.6(-ll) 


C 9 H 




4.1(-14) 


7.2(-ll) 


2.9( 


-11) 


2.3(-ll) 


HC 5 N 


9.7(-ll) 


1.4(-12) 


1.6(-10) 


1.4( 


-10) 


1.5(-10) 


HC 7 N 


2.7(-ll) 


7.5(-14) 


2.9(-ll) 


2.1' 


-11) 


2.0(-ll) 


HCgN 


2.5(-12) 


5.0(-15) 


5.9(-12) 


3.0 


-12) 


2.7(-12) 


C 2 H 2 




6.5(-09) 


3.7(-08) 


3.3( 


-08) 


3.2(-08) 


C 3 H 2 


l.O(-ll) 


5.0(-13) 


l.l(-ll) 


9.6( 


-12) 


9.3(-12) 


C4H2 


2.2(-ll) 


2.4(-12) 


1.2(-10) 


9.6( 


-11) 


9.0(-ll) 


C5H2 




1.3(-13) 


2.5(-ll) 


1.8( 


-11) 


1.6(-11) 


C 6 H 2 


2.5(-12) 


6.3(-14) 


2.5(-ll) 


1.5( 


-11) 


1.4(-11) 


C 7 H 2 




9.8(-15) 


7.1(-12) 


4.0( 


-12) 


3.4(-12) 


C 8 H 2 




7.1(-15) 


5.6(-12) 


2.6( 


-12) 


2.1(-12) 


CgH 2 




5.5(-16) 


2.6(-12) 


1.0( 


-12) 


7.8(-13) 


CH3CCH 


7.8(-10) 


5.9(-12) 


2.2(-09) 


2.1( 


-09) 


2.2(-09) 


CH3C4H 




1.7(-12) 


1.1 (-09) 


7.0 


-10) 


6.2(-10) 


CH 3 C 6 H 




3.1(-13) 


1.7(-10) 


8.6( 


-11) 


6.8(-ll) 


CCS 


8.5(-ll) 


8.2(-ll) 


1.3(-10) 


1.4' 


-10) 


1.4(-10) 


HCO^ 


1.0(-12) 


1.4(-12) 


4.1(-13) 


3.8( 


-13) 


3.7(-13) 


HCN 


1.2(-09) 


4.6(-10) 


7.7(-10) 


8.3( 


-10) 


9.3(-10) 


HNC 


3.2(-10) 


3.9(-10) 


6.3(-10) 


7.0( 


-10) 


7.8(-10) 


HC 3 N 


8.9(-10) 


2.0(-ll) 


4.0(-10) 


4.3( 


-10) 


4.7(-10) 


CO 


3.9(-05) 


8.9(-06) 


5.1(-05) 


5.4( 


-05) 


5.4(-05) 


CS 


3.3(-10) 


6.4(-10) 


5.1 (-09) 


5.1 


'-09) 


4.6(-09) 


SO 


1.4(-10) 


9.0(-10) 


1.3(-10) 


2.1( 


-10) 


1.9(-10) 


HCO+ 


6.0(-10) 


4.1(-10) 


1.7(-09) 


1.7' 


-09) 


1.8(-09) 


N 2 H+ 


2.5(-10) 


2.7(-ll) 


l.l(-ll) 


l.Oi 


-11) 


l.O(-ll) 


CN 


8.0(-ll) 


1.4(-10) 


9.6(-10) 


1.2 


'-09) 


1.3(-09) 


C 2 




2.0(-10) 


l.l(-09) 


1.1( 


-09) 


l.l(-09) 


C 3 




1.5(-07) 


l.l(-08) 


7.4( 


-09) 


6.1(-09) 


C4 




2.8(-ll) 


1.0(-09) 


8.2( 


-10) 


6.9(-10) 


C 8 




1.3(-08) 


6.7(-09) 


3.7( 


-09) 


2.7(-09) 


C e 




2.1(-11) 


l.l(-09) 


5.8( 


-10) 


4.4(-10) 


C7 




6.0(-10) 


2.7(-09) 


1.1 


-09) 


7.4(-10) 


C 8 




9.7(-13) 


4.3(-10) 


1.5 ( 


-10) 


1.1(-10) 


C 9 




3.4(-ll) 


1.3(-09) 


3.5( 


-10) 


2.1(-10) 


C10 




3.2(-15) 


7.3(-ll) 


1.8( 


-11) 


1.1O-11) 


C 2 N 




9.9(-ll) 


4.3(-10) 


4.0( 


-10) 


3.9(-10) 


C 3 N 




1.9(-11) 


1.2(-10) 


1.2' 


-10) 


1.4(-10) 


C 4 N 




9.9(-12) 


1.7(-10) 


8.2( 


-11) 


6.3(-ll) 
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Tabic 5 — Continued 





Tmax - 


= 10 K 


^max - 


= 30 K 


Tmax — 


100 K 


Tmax - 


= 200 K 


Species 


Obs. a Xi 


1:* p t. 




1:* p t. 


x t 


O p t 


Xi 


,opt 


C 5 N 


5.9 


-13) 


1.1 


-10) 


1.2 


-10) 


\.i 


-10) 


C 7 N 


3.8 


-14) 


2.6 


-11) 


1.6 


-11) 


1.5 


-11) 


C 9 N 


2.2 


-15) 


5.1 


-12) 


1.2 


-12) 


1.9 


-12) 


CH3C3N 


1.9 


-14) 


2.9 


-13) 


2.7 


-13) 


2.7 


-13) 


CH3C5N 


1.7 


-16) 


2.3 


-14) 


1.9 


-14) 


1.8 


-14) 


CH3C7N 


8.0 


-18) 


3.8 


-15) 


2.3( 


-15) 


2.2 


-15) 


N 2 


7.6 


-07) 


5.6 


-06) 


5.7 


-06) 


5.7 


-06) 


O2 


3.6 


-07) 


1.1 


-08) 


5.1 


-09) 


4.6 


-09) 


CO2 


1.7 


-07) 


4.1 


-09) 


2.8( 


-09) 


2.4 


-09) 


C3H3 


2.3 


-12) 


2.2 


-09) 


2.2( 


-09) 


2.3 


-09) 


C4H4 


4.1 


-12) 


9.5 


-11) 


7.8 


-11) 


8.5 


-11) 


C 6 H 6 


3.0 


-13) 


3.2 


-10) 


1.8 


-10) 


1.7 


-10) 


H 2 CS 


1.4 


-11) 


5.2 


-10) 


6.6 


-10) 


7.3 


-10) 


H 2 S 


1.1 


-10) 


3.3 


-09) 


7.8 


-09) 


7.7 


-09) 


CH 2 CN 


2.7 


-11) 


7.4 


-11) 


8.6( 


-11) 


8.6 


-11) 


CH2CO 


2.7 


-10) 


1.0 


-09) 


1.2 


-09) 


1.1 


-09) 


CH3CN 


6.3 


-13) 


2.0 


-12) 


2.3( 


-12) 


2.6 


-12) 


CH3OH 


2.9 


-12) 


5.0 


-11) 


3.4 


-12) 


1.8 


-12) 


CH3OCH3 


1. 1 


-17) 


3.8 


-16) 


1.1 


-17) 


6.9 


-18) 


H 2 CO 


5.9 


-11) 


4.3 


-09) 


4.3 


-09) 


4.5 


-09) 


HCOOH 


6.1 


-11) 


3.6 


-12) 


1.7 


-12) 


1.5 


-12) 


HCOOCH3 


2. 1 


-16) 


3.8 


-15) 


5.9( 


-16) 


4.7 


-16) 


HNCO 


2.2 


-10) 


2.4 


-12) 


4.3( 


-12) 


4.2 


-12) 



a Obscrvational references are the same as those in Table [3] 

b Boldface type indicates that the computed abundance is less than the observed value by 
one order of magnitude or more. 

c Italic type indicates that the computed abundance is greater than the observed value by 
one order of magnitude or more. 



